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We propose a new description of astronomical dust emis- 
sion in the spectral region from the Far-Infrared to millimeter 
wavelengths. Unlike previous classical models, this descrip- 
tion explicitly incorporates the effect of the disordered inter- 
nal structure of amorphous dust grains. The model is based 
on results from solid state physics, used to interpret laboratory 
data. The model takes into account the effect of absorption by 
Disordered Charge Distribution, as well as the effect of absorp- 
tion by localized Two Level Systems. We review constraints 
on the various free parameters of the model from theory and 
laboratory experimental data. We show that, for realistic values 
of the free parameters, the shape of the emission spectrum will 
exhibit very broad structures which shape will change with the 
temperature of dust grains in a non trivial way. The spectral 
shape also depends upon the parameters describing the inter- 
nal structure of the grains. This opens new perspectives as to 
identifying the nature of astronomical dust from the observed 
shape of the FIR/mm emission spectrum. A companion paper 
will provide an explicit comparison of the model with astro- 
nomical data. 



1. Introduction 

It is now well established that the Spectral Energy Distribution 
(SED) from the Inter-Stellar Medium (ISM) emission from our 
Galaxy and most external galaxies is dominated by thermal 
emission from dust grains which spans over almost 3 orders 
of magnitude in wavelengths, from the Near Infra-Red (NIR) 
to the millimeter wavelength range. Several dust models have 
been proposed to explain the observations. Recent dust models 
share many common characteristics (e.g. Mathis at al. 1 1 9771 
Draine & Lee 119841 Draine & Anderson Q 9851 Weiland et al. 
[19861 Desert et al. [19901 Li & Greenberg [1997] Dwek et al. 
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119971 Draine & Li 2001 ). Most of them require a wide range 
of the dust grain sizes spanning almost 3 orders of magnitude, 
from sub-nanometer dimensions to fractions of a micron. They 
also require a minimum of three components with different 
chemical composition in order to explain the ultraviolet and 
optical extinction, along with infrared emission. 

The smallest component is needed to explain the "aro- 
matic" emission features at 3.3, 6.2, 7.7, 8.6 and 11.3yum. 
It is now routinely assigned to large aromatic molecules or 
Polycyclic Aromatic Hydrocarbons (PAH) (Leger & Puget 
[19M Allamandola et al. [19831 Desert et al. [19901 PAH are 
transiently heated after absorbing single UV and far-UV pho- 
tons, and cool down emitting NIR photons in vibrational tran- 
sitions representative of their aromatic structure. 

A second component, composed of Very Small Grains 
(VSG) is required to explain the continuum emission in the Mid 
Infra-Red (MIR). Emission in this range still requires signifi- 
cant temperature fluctuations of the grains, which necessitates, 
under conditions prevailing in the Inter-Stellar Medium (ISM), 
grains sizes in the nanometer range. VSGs could be composed 
of carbonaceous material whose absorption could also explain 
the 2200A UV bump in the extinction curve (see for instance 
Desert et al. [l990V 

A third component, composed of Big Grains (BG), with 
sizes from 10-20nm to about lOOnm, is necessary to account 
for the long wavelength emission (in particular as observed in 
the IRAS 100 jum band and above). The BG component domi- 
nates the total dust mass and the absorption in the visible and 
the NIR. Observations of the "silicate" absorption feature near 
10 /urn indicate that their mass is dominated by amorphous sil- 
icates (see Kemper et al. 120041 1 . Such amorphous structure for 
these dust grains is expected from the interaction with cosmic 
rays which should alter the nature of the solid (Brucato et al. 
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2004, Jager et al. 2003 ), but also from the formation processes 
of these BGs by aggregation of sub-sized particles. 

Recent studies attempting to explain the submillimeter ob- 
servations, revealed the need for either the existence of an ad- 
ditional most massive component of very cold dust (e.g. Reach 
et al. fl995l Finkbeiner et al. [19991 Galliano et al. 120031 or for 
substantial modifications of dust properties in the submillime- 
ter with respect to predictions of standard dust models (e.g. 
Ristorcelli et al. [T998l Bernard et al. [l999l Stepniket al. l2003al 
Dupac et al. 120011 Dupac et al. 120021 Dupac et al. 120031 . 

Similarly, laboratory spectroscopic measurements of inter- 
stellar grains analogs reveal that noticeable variations in their 
optical properties in the Far Infra-red (FIR) and Millimeter 
(mm) wavelength range can occur (Agladze et al. 1994 
Agladze et al. [19961 Mennella et al. [I99"8l Boudet et al. [2005l . 
They pointed out the possible role of two phonons differ- 
ence processes, disorder-induced one-phonon processes, res- 
onant or relaxation processes in the presence of Two-Level 
Systems(TLS) in amorphous dust (see Sect. [43). The influence 
of such Two-Level Systems on interstellar grain absorption and 
emission properties was first proposed to the astronomers's 
community by Agladze et al. 119941 More recently, a prelim- 
inary investigations by Boudet et al. 2002 showed that opti- 
cal resonant and relaxation transitions in a distribution of TLS 
could qualitatively explain the anticorrelation between the tem- 
perature and the spectral index of dust as observed by the 
PRONAOS balloon experiment (e.g. Dupac et al. 120031 . 

A precise modeling of the long wavelength dust emission 
is important in order to accurately subtract foreground emis- 
sion in cosmological background anisotropy measurements, 
especially for future missions for space CMB measurements 
(Lamarre et al. 120031) and for surveys of foreground compact 
sources (see Gromov et al. 12002). A number of comprehensive 
foreground analysis papers (Bouchet & Gispert 1999 Tegmark 
et al. 120001 Bennett et al. l20"03l Barreiro et al. 120061 Naselsky 
et al. 2005 and other) did not pay attention to dust emissivity 
model accuracy. In addition, precise modeling of the FIR/mm 
dust emission is also important to derive reliable estimates of 
the dust mass, to trace the structure and density of pre-stellar 
cold cores in molecular clouds. Indeed the dust governs the 
cloud opacity, and influences the size of the cloud fragments. 
In addition, the efficiency of star formation should be related to 
the dust properties, which are expected to vary between the dif- 
fuse medium and the denser cold cores. It is thus of importance 
to know how some physical or chemical properties of dust can 
influence their FIR/mm emission. 

For models of the submillimeter electromagnetic emission 
of dust, considering the intrinsic mechanical vibrations of the 
grain structure is natural since they fall in the right frequency 
range (typically in the frequency range v ~ 10 11 - 10 13 Hz, i.e. 
A ~ 30 jum - 3 mm). Observations in this spectral region poten- 
tially provides new tools for investigating the internal structure 
of dust grain material. The internal mechanics of grains is also 
important for chemistry in the ISM. The dust vibrations sup- 
ply an energy sink for newly formed molecules which are usu- 
ally formed in unstable excited states. Understanding the grain 
structure and intrinsic movements is necessary for considera- 



tions regarding grain collisions, destruction and agglomeration, 
including preplanetary bodies formation in circumstellar disks. 

In Sect. [21 we first recall some basics regarding the FIR/mm 
dust emission and the semi-classical model of light interac- 
tion with matter. In Sect. [3) we gather evidences for spectral 
variations in astronomical and laboratory measurements. Then, 
we present in Sect. [4] a model, based on the intrinsic proper- 
ties of amorphous materials, taking into account absorption by 
Disordered Charge Distribution as well as the effect of absorp- 
tion by a distribution of TLS. In Sect.[5] we propose a new 
model of the submillimeter and millimeter grain absorption and 
discuss the range of plausible values for the free parameters, 
based on existing laboratory data. In Sect. 15.41 we discuss the 
implications of the model, in particular regarding the expected 
variations of the FIR/mm emission spectra with the dust tem- 
perature. Sect. [6] is devoted to conclusions. The determination 
of the model parameters applicable to astronomical data will 
addressed in a companion paper. 

2. Basic knowledge on the FIR/mm dust emission 

2.1. Interstellar dust emission and extinction 

The intensity of thermal emission from interstellar dust at tem- 
perature 7d is 

I b) {io,T) = e e {co)-B <a {io,T i ) (1) 

where 1^ is energy flux density per unit area, angular frequency 
and solid angle (or specific intensity), e e the dust emissivity and 
B w the Planck function at angular frequency to. 

According to the Kirchhoff law, the emissivity e e is equal 
to the absorptivity 

e e (cj) = l-e- T ^\ (2) 

where the optical depth r is related to the dust mass column 
density on the line of sight M c { and the dust mass opacity k as 

t(oj) = k((jS) ■ M d (3) 

For spherical grains of radius a and density p, the dust opac- 
ity (effective area per mass) is given by 

k(u) = , (4) 

4p a 

where the absorption efficiency Q(oS) = cr(cu) / (na 2 ) is the ratio 
of the absorption cross section cr(a>) to the geometrical cross 
section of the grain na 2 . The grain equilibrium temperature 
is deduced from the balance between the emitted and absorbed 
radiation from the Inter-Stellar Radiation Field (ISRF) 

Jr»oo /-»oo 
BJco, T d )Q(Lo)dco = ll SRF Q(aj)dcj (5) 
o Jo 

The Mie theory, using the Maxwell equations of the elec- 
tromagnetic theory, leads to an exact solution for the absorption 
and scattering processes by an homogeneous spherical particle 
of radius a whose material is characterized by its complex di- 
electric constant e(a>) = e'(cj) + ie"(ai). In that case, the ab- 
sorption coefficient Q(oj) writes as an infinite series of terms 
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related to the complex dielectric constant. For particles small 
compared with the wavelength such as | Ve(w)| • 2na/A <K 1 
(where A = 2nc/o)is the wavelength and c the speed of light in 
vacuum), the absorption coefficient can be approximated by 



Q(co) = — • Im 



e(a>) - 1 



e(w) + 2 



(6) 



In the FIR/mm wavelength range, for nonconducting mate- 
rials of astronomical dust particles , the condition e" <K e' is 
generally satisfied, and in this case the absorption coefficient 
reduces to 
Q(to) oj 12 



(7) 



a c (e' + 2) 2 ' 
As a consequence, following Eq. (0 the dust opacity 



rewrites 



k(oj) = 



c (e' + 2) 2 p ' 



(8) 



2.2. Assumption of a A-independent emissivity 
spectral index 

In the absence of specific knowledge about realistic dependen- 
cies for e(oj), a simple power law approximation for FIR/mm 
dust emission is often assumed 



co 



(9) 



where j3 is referred to as the emissivity spectral index. 

The intensity of thermal emission from interstellar dust, 
assuming an optically thin medium, at a given wavelength 
A = 2jt/(l> located in the FIR/mm range is 



I u (<o, T) = £ e (Ao) ■ — 
Wo 



BJu,T d ) 



(10) 



where e e (Ao) is the dust emissivity at a reference wavelength 
Aq — 2ttc/loo. 

Simple semi-classical models of absorption such as the 
Lorentz model for damping oscillators and the Drude model 
for free charge carriers (Sect. l2.4i i provide an asymptotic value 
P - 2 when co — > 0. Such a value of the spectral index was 
in satisfactory agreement with the earliest observations of the 
FIR/mm SED of interstellar dust emission, but not necessarily 
with the most recent ones. 

The emissivity spectral index j3 is equal to the slope of a 
plot of the dust opacity versus wavelength in logarithmic scale. 
In the general case, the power-law assumption is not valid, and 
the spectral index has to be considered as a wavelength and 
temperature dependent parameter. Similarly, e e (^o) in Eq. ( TTOb 
should be considered a function of T. 

2.3. Relations between absorption properties of bulk 
material and dust 

The absorption coefficient defined as the optical depth in the 
bulk material per unit length is 



a = 2lm(k), 



(11) 



where k is the complex value of the wave number k — a> -\fe/c 
and e is the complex dielectric constant Therefore for e" <K e', 



(12) 



In dust and porous material measurements, a mass absorp- 
tion coefficient a /pis also often used. 



In the following analysis (Sect. 12.41 and other) we will link 
the macroscopic value e to the value of the local susceptibility 
X - x' + i-X'o defined on microscopic scale. 

Because small astronomical particles cannot always be con- 
sidered as a continuous medium, the electrodynamic equations 
of continuous media must be applied with care. So we detail in 
the following the equations used, and we discuss their range of 
application. 

The Equations of motion describes the response of charges 
to a local electric field Eq, which permits to determine a dipole 
moment per unit volume as P — xqEo- It should be taken into 
account that the local electric field Eq differs from an external 
field E due to additional field produced by other parts of the 
dielectric. In the case of isotropic dielectric materials, 

E = E + —P. 

The dielectric constant e links the external field E and the elec- 
tric inductance D — E + AnP through D = eE. In this case, 
e can be derived from the susceptibility Xo, using the so called 
Clausius-Mossotti equation 



e = 1 + 



1 - f^o 



When e" <K e', we have 



(e'+2) 2 



Wo, 



(13) 



(14) 



which leads to 



Anu (e' + 2) 2 

a(aj) = — — ^(w, 

cVe 7 9 



k(oj) — 



c p 



(15) 
(16) 



In Eq. ([T5ll. the term (e' + 2) 2 /9 is the local field correction 
factor, the term Ve 7 is the refraction index. 

Thus, the relation between the absorption coefficient a de- 
duced from transmission measurements in the laboratory and 
the absorption efficiency Q or opacity k, for small spherical 
grains of radius a and density p can be obtained from equations 
Eq.l7land[T2l and writes 



Q(u>) _ 12a- VF 
a(u) ~ (e' + 2) 2 

and as a consequence 

k(oj) _ 9 ■ yf? 
a(a>) ~ p ■ {P + 2) 2 



(17) 



(18) 
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The condition e" <k e' is valid for dust materials in the 
FIR/rnm wavelength range, implying that the real part of the 
dielectric constant e' can be considered to first order as a con- 
stant, leading to de' /da> = 0. This comes from the coupling 
between the real part e' and the imaginary part e" through 
the Kramers-Kronig relations exposed in Sect. 13.21 This implies 
that, in that spectral range, the absorption coefficient a, the ab- 
sorption efficiency Q and the opacity k are simply proportional. 
Thus the slopes of a(co), Q{a>) and k{u) are the same. As for the 
emissivity, an absorptivity spectral index can then be defined by 
the slope of the spectral plot of the absorption a(u>) in logarith- 
mic scale. This absorptivity spectral index deduced from labo- 
ratory transmission measurements can be directly compared to 
the emissivity spectral index deduced from astrophysical mea- 
surements. 

2.4. The semi-classical models 

The well known semi-classical model of light interaction with 
charge carriers in matter (electrons and ions) considers the mo- 
tion of charged particles driven by the time-dependent electric 
field of a monochromatic plane wave, magnetic forces being 
neglected compared with electrical forces. 

The equations of motion of interacting charge carriers are: 



(19) 



re 



where mj, qj are masses and charges of particles, Kjt are con- 
stants of elastic interaction, is ^-th component of electrical 
field vector E, w. are components of the 3 -dimensional vec- 
tor of /'-th ion displacement. Damping parameters bj describes 
dissipation of energy. Force constants K are the coefficients of 
power expansion of the potential energy U. 



U = l V Kf'i/J., + ... 

9 Z_i J J J J 

re 



(20) 



U(- ■ ■ Uj ■ • ■) depends only on relative deviations (u - u\) 
with corresponding symmetry of Ky., , particularly 



A solution of Eq. (19$ with E(t) = £}e~'°" is u = ue 
qj E ■ e* Kj u>l 



K J 



u 2 - a) 2 - iy K (x> 



(24) 



Dielectric properties of a substance can be described by a 
susceptibility xo, which is the dipole moment per unit volume 
generated by a unity local strenght = 1 electrical field, 



p K u K 



o>£ - or - iy K to 
For isotropic materials, constants 



(25) 



(26) 



does not depend upon the direction of E. 

In the extreme case where oj k » y K the spectrum of the 
dielectric absorption can be described as a collection of non- 
overlapping resonant Lorentz profiles 



Xo = 



(«2 - w 2 ) 2 + (y K to) 2 ' 



(27) 



It is the case of infrared absorption bands, especially in 
ionic crystals ("reststralen" bands). This corresponds to the 
third term of Eq. (fT~9b being small. 

In the opposite case {lo k <k y K ) this leads to a Debye relax- 
ation spectrum when co <k w ( : 



„ _ p K (ojt k ) 

X ° ~ 1+(WT,) 2 ' 



(28) 



where t k - y K lco 2 K is a relaxation constant. This is usually the 
case of polar dielectric liquids. It corresponds to the first term 
of Eq. ( fT9l being negligible. 

This spectrum corresponds to situations where the inertial 
effects (first term in Eq.[T9l) are comparatively small. If this 
condition is satisfied, the time constant in Eq. d28T i can usu- 
ally be determined by measuring or calculating the timescale 
of energy dissipation during return to equilibrium state, with- 
out solving Eq. ( fT9l . 

A third extreme case, corresponding to the second term of 
Eq. (f28b being negligible is usually refered to as the Drude 
approximation. It applies to free carriers and is considered in 
Sect.l2~6l 



K 



7'f 



(21) 2.5. The Lorentz model for crystalline insulators 



Linear Eq. (fT9] > with E — have solutions in the form of a 
sum of damping oscillations 

where 

« = (••• M w • • •) is the multidimensional deviation vector, a K 

are constants, e K = (• • • e\ ■ ■ •) are eigenvectors, and y K -iu K are 
eigenvalues of matrix 



(23) 



The Lorentz model applies for crystalline isolators with an 
ionic character, and describes the resonant interaction of an 
electromagnetic wave with the vibrations of ions which act like 
bound charges. Such an interaction is caused by the electric 
dipoles which arise when positive and negative ions move in 
opposite direction in their local electric field. This kind of dust 
absorption can be clearly identified in astrophysical observa- 
tions as emission or absorption bands. Indeed, the optical depth 
t is proportional to opacity k (Eq.fToT): 



k(oS) 



_ 47r v 



p K a? K y K o? 



u 2 ) 2 + (y K u) 2 



(29) 
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The maximum number of optically active vibrational 
modes in a crystal is n opt - 3(s - 1), where s is the number 
of atoms in an elementary cell of the lattice. They are usually 
classified as branches in dispersion dependence of oscillation 
frequency u)(k) on wave number k. All these branches are re- 
stricted by relatively high a>(k) in the infrared region. Another 3 
branches with a>(k) — > when k — >0 are acoustic modes. They 
are optically not active in crystal because a total charge of cell 
is zero. Another situation for amorphous matter is considered 
here in Sect.l4~2l 

For substances of astrophysical interest, the characteristic 
frequencies of bands fall in region A £ 50 pm. Therefore, in the 
FIR/mm wavelength range, the asymptotic a> — > of Eq. d29l ) 
can be used: 

k{u>) = u L — ) —r- (30) 

Thus, the Lorentz model formally applied to dust in the 
FIR/mm wavelength range gives a spectral index jS = 2, a value 
often used in "standard" dust models. 

Many dust models have related the long wavelength absorp- 
tion by interstellar silicates and ices (A > 100 pm) to the damp- 
ing wing of the infrared-active fundamental vibrational bands, 
visible at shorter waves (see for example Andriesse et al. 1 19741 
and references therein). 

Laboratory infrared data are satisfactorily described (see 
for example Rast 1968) by a sum of resonant Lorentz profiles 
(or reststrahlen bands) as described by Eq. d25l l and Eq. ( f29l . 
in which the various constants p K , a> K , y K are determined from 
spectroscopic measurements. This model gives a satisfactory 
approximation for the infrared region, not only for ionic crys- 
tals, but also for covalent ones and for their amorphous coun- 
terparts with an appropriate choice of p K , a> K , y K . 

The situation is more complicated for the low frequency 
region, where absorption is defined by damping constant y, in 
contrast to IR bands, where absorption is weakly dependent on 
y. The damping constants bj introduced in the semi-classical 
model (Eq.[T9ll and therefore y K in Eq. ( 1291 have no simple 
microscopic explanation. The damping is linked to heating or 
distribution of energy to all modes of grain vibrations. So y 
characterizes intermode interactions which are frequency and 
temperature dependent in general. In contrast to single phonon 
absorption of resonant IR photons, a wing absorption is a two 
phonon difference process. Due to the conservation of energy, 
the energy of the absorbed low frequency photon is equal to 
the energy difference of the two phonons in the fundamen- 
tal vibrational bands (Rubens and Hertz 119121 Ewald 119221 . 
In fact in crystals, any combination of acoustical and optical 
phonons which satisfies the energy and momentum conserva- 
tion rules can take part to the absorption . This is the main 
source of infrared absorption in homoatomic crystals. In the 
FIR wavelength range the two-phonon difference processes are 
responsible for the absorption in ionic crystals. Because the ef- 
ficiency of such processes are entirely determined by the dif- 
ference between the phonon occupation numbers of the two in- 
volved phonons, the corresponding absorption decreases dras- 
tically at decreasing temperatures (Hadni [T9701 l. 
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However the optical behavior of amorphous materials is 
different. The disorder of the atomic arrangement leads to a 
breakdown of the selection rules for the wavenumber. The 
absorption spectra reflects thus the whole density of vibra- 
tional states. This induces in the longest wavelength range a 
broad absorption band due to single phonon processes which 
completely dominates over the multiphonon absorption (H. 
Henning and H. Mutschke 1997). As a consequence, such lat- 
tice absorption does not vanish at low temperature for the amor- 
phous solids , in contrast to the crystalline ones. 

2.6. The Drude model for conducting materials 

For free charge carriers the previous semi-classical equa- 
tion of motion (Eq.[T9l can still be used, setting all restor- 
ing force coefficients K to zero, equaling all the masses 
rrij to the effective mass of the free carriers m c ( mj 
m, ). and using a constant value for the damping parameters 
bj/mj=const=y c . Therefore Eq. (f25j) for the susceptibility sim- 
plifies as Xo — - Oi>p (cl> 2 - iy c (Sf X < where u> p is the so-called 
plasma frequency to 2 = N c q 2 c lm c , and N c is the number of free 
carriers per unit volume. 

In metals, the number of free electrons per unit volume is 
independent of temperature, and is of the same order of mag- 
nitude than the number of atoms (N c ~ 10 22 cm 4 ). This leads 
to a frequency o> p ~ 10 15 s located in ultraviolet or in the 
visible, and to correspondingly high conductivity cr = to 2 p ly c . 

In the long wavelength range, far from the plasma fre- 
quency (<d <K (D p ) and when the modulus of e is large compare 
to 1, the asymptotic behavior of the opacity is given by 

k(w) = ^ 2 -^L (31) 
cpof p 

which also gives a spectral index /3 = 2. Metallic materials and 
high conductivity semiconductors are expected to follow this 
model. Possible dust absorbents of such type are graphite and 
metallic iron. 

Therefore, both the Lorentz and the Drude models in the 
low frequency limit yield ji-2 in Eq. (0. These very different 
models have common features. Their characteristic times (u>~ 1 , 
o)p l ) are far from FIR/mm region. In this case a good approx- 
imation is obtained by asymptotic expansion of the absorption 
coefficient a{a)) on even powers of frequency oj 2 ", n > 0. In the 
general case the first term is dominating (n — 1) and then we 
have P-2. 

Other values (e.g. n-2, y8=4) are possible if physical mech- 
anisms suppressing the n — 1 term are present, as in the case 
of screening considered in Sect. 14. 21 Intermediate values of ji 
provide a clear evidence for specific processes which time con- 
stants are of the order of aT 1 . 

3. Temperature and spectral variations of dust 
optical properties 

3.1. Astronomical evidences 

Analyses of millimeter and submillimeter emission from 
molecular clouds have found spectral indices between /J=1.5 
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(Walker, Adams & Lada,1990) and /3=2 (Gordon 1988,Wright 
et al., 1992). However some values in excess of 2 can also be 
found in the literature (Schwartz 1982, Gordon 1990). A few 
studies of dense clouds have yielded spectral indices around 
1 (Gordon 1988, Woody et al. 1989, Oldham et al. 1994), but 
observations of disks of gas and dust around young stars can 
indicate values sometimes less than 1 (Chandler et al. 1995). 
A submillimeter continuum study of the Oph 1 cloud (Andre, 
Ward-Thompson & Barsony 1993) found a typical value of /3 
equal to 1 .5 but attributed changes in observed ratios to tem- 
perature variations. Without appropriate temperature data, the 
authors could not be conclusive on this issue. 

More recent observations have evidenced that the actual 
spectral energy distribution of dust emission could be signifi- 
cantly more complicated than described above. Analysis of the 
FIRAS results have shown that, along the galactic plane, the 
emission spectrum is significantly flatter than expected, with a 
slope roughly compatible with /3 = 1 (see Reach et al. 1995 
in particular their Fig. 7). They attributed the flattening com- 
pared to the /3 — 2 canonical value to an additional component 
peaking in the millimeter range and favored that the millime- 
ter excess could be due to the existence of very cold dust at 
Ti = 5 - 7 K in our Galaxy. However, the origin of such very 
cold dust remained unexplained, and the poor angular resolu- 
tion of the FIRAS data (^ 7°) did not allow further investiga- 
tions. 

Finkbeiner et al. [1999 later adopted a similar model. They 
showed that a good fit of the overall FIRAS data could be ob- 
tained using a mixture of warm dust at 7^ =s 16 K and very cold 
dust at 7^ =! 9 K. The temperature of the warm component was 
set using a fit to the DIRBE data near A = 200 fim while the 
temperature of the very-cold component was set assuming an 
independent dust component immersed in the same radiation 
field but with different absorption properties. They implied that 
these two components could be graphite (warm) and silicates 
(very cold). However, they gave no quantitative argument to 
support this hypothesis 

PRONAOS was a french balloon-borne experiment de- 
signed to determine both the FIR dust emissivity spectral in- 
dex and temperature, by measuring the dust emission in four 
broad spectral bands centered at 200, 260, 360 and 580 fim (see 
Laman-e et al ll994l Serra at al 120021 Pajot et al l2006t . The 
analysis of the PRONAOS observations towards several re- 
gions of the sky ranging from diffuse molecular clouds to star 
forming regions such as Orion or M17 have evidenced an anti- 
correlation between the dust temperature and the emissivity 
spectral index (see Dupac et al. 2003). The dust emissivity 
spectral index varies smoothly from high values (ft =s 2.5) for 
cold dust at Td - 12 K to much lower values (J3 — 1) for dust 
at higher temperature (up to 80 K) in star forming regions 
like Orion. Dupac et al. ( 120031 ) showed that these variations 
were not caused by the fit procedure used and were unlikely to 
be due to the presence of very cold dust. They concluded that 
these variations may be intrinsic properties of the dust. 

Extragalactic observations have also revealed unexpected 
behavior of dust emission at long wavelength. For instance, 
Galliano et al. ( l2003ll2003T l. combining JCMT data at 450 and 
850 pm, IRAM 30m data at 1.2 mm with ISO and IRAS mea- 



surements in the IR and FIR evidenced a strong millimeter ex- 
cess towards a set of four blue compact galaxy (NGC1569, II 
Zw40, He 2-10 and NGC 1 140). They attributed this excess to 
the presence of very cold dust (7^ - 5 -1 K) with a very flat 
emissivity index (J3 — 1) which would have to be concentrat- 
ing into very dense clumps spread allover the galaxy. This very 
cold dust would then account for about half or more of the total 
dust mass of those galaxies. 

3.2. Laboratory evidences 

3.2.1. Comparing laboratory and astronomical data 



As it was shown in Sect. |2.3l laboratory data on absorption 
spectral index could be used for interpreting astronomical ob- 
servations, taking into account the quantitative limitations con- 
sidered here. Some other limitations should be considered. 

First, there is a difference between bulk laboratory samples 
and the small cosmic dust. The latter probably has a high sur- 
face to volume ratio, and surface effects are important for spec- 
troscopy as for the physics and the chemistry of dust. 

Second, samples used for laboratory measurements are 
generally synthesized in very small amounts. Along with the 
small values of the absorption coefficients it requires high- 
sensitivity dedicated instruments such as 3 He-cooled bolome- 
ters or/and high power sources of radiation. 

In the laboratory indirect IR/FIR/mm spectroscopy meth- 
ods are also used, such as Kramers-Kronig spectroscopy. The 
dielectric constant counterparts (e' - 1) and e" are Hilbert 
transforms of each other and are related through the Kramers- 
Kronig relations (see Landau & Lif shits [1 9821 > by 



e'ipS) 



e"(x) 



i r° 

n J-o. 



= _i r eX 

n J-co x - 



dx + 1, 



(32) 



dx, 



where both integrals are the Cauchy principal values. 
A simple form of Eq. (l32l is for o> — 0: 

2 r° „ 

e'(0) = - e"(x)dlnx+ 1. 
t Jo 



(33) 



The use of these relations together with reflection measure- 
ments R(u>) = |(V^ _ 1)/(V^ + l)l 2 allows to resolve the sys- 
tem of equation with respect to e'(w) and e"{o>) and therefore 
to calculate an absorption spectrum. Possible method inaccu- 
racies due to inverse problem solution should be considered 
carefully. In particular, it is sometimes assumed as a prior that 
the solution has the form of a power spectrum. In practice, the 
most reliable approach should be to use the same model for 
interpretation of both the astrophysical and laboratory data. 

3.2.2. Spectral index variations in disordered solids 

It has been known for more than 25 years, that disordered solids 
(glasses and mixed crystals) have significantly higher submil- 
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limeter absorption than perfect crystals of similar chemical na- 
ture. 

Mon et al. ( 1 19751 ) performed first absorption measurements 
on various amorphous dielectric materials in the millimeter 
wavelength range (1 mm - 5 mm) and at low temperature (in 
the range 0.5 K - 10 K). A strong temperature-dependence of 
the absorption was observed, leading to an millimeter absorp- 
tion excess at low temperature. They already attributed this ef- 
fect to the presence of a distribution of Two-level Systems in 
the studied materials. 

Others experimental results on bulk materials at frequen- 
cies to line between 0.1 and 150 cm -1 and T between 10 K and 
300 K were discussed by Strom & Taylor (1977I ). The spectra 
follow a general empirical relation a = Ka> 2 between 15 and 
150 cm 1 , where the temperature independent constant K de- 
pends on the material considered and its internal structure. At 
lower frequencies w/27r <15cm _1 (A > 700 /mi) the slope of 
the spectra is variable (J3 between and 3) and temperature de- 
pendent as shown in Strom & Taylor ( 11 9771 1 . 

Bosch ( 1978 ) performed absorption measurements for tem- 
perature between 300 K and 1.6K in the FIR/mm range 
(500 /mi - 10 mm) on amorphous silicates mainly composed 
of SiC>2. The studied silicate is not a typical grain analog. 
However, measurements performed over such a wide spectral 
and temperature range are particularly appropriated to study 
physics governing absorption processes. Bosch found a strong 
temperature and frequency dependency of the absorption coef- 
ficient in the millimeter range (A > 500 /mi). The absorption 
coefficient was characterized by an increase of the spectral in- 
dex with temperature from /3 = 1.6 at 300 K to y8 = 3 at 10 K 
and by a strong temperature dependence. To describe this tem- 
perature and frequency behavior, Bosch also referred to the ex- 
istence of Two Level Systems in the material. 

Therefore, laboratory spectroscopy shows evidences for de- 
viation of absorption spectral index from the canonical f3 — 2 in 
the FIR/mm. The temperature independent absorption process 
which provides ft = 2 in the spectral region A < 700 /im does 
not seem to dominate the absorption at longer wavelengths. We 
will see in Sect. |4] that the model proposed here includes this 
phenomenon in a natural manner. 

3.2.3. Spectral variations in dust analogues 

FIR/mm optical properties of solids with chemical composition 
and structure considered as analogues of interstellar dust were 
investigated in a few laboratory studies (for example Koike 
et al. [19801 Agladz e [19961 , He nning and Mutschke [19971 
Mennella et al. 119981 Boudet et al. 12005) . Particular interest has 
been given to low temperature investigations of amorphous or 
other disordered materials, which reveal systematic deviations 
from the phonon theory of crystal solids and was interpreted 
using the TLS theory (see Bosch [l978l Phillips [1987] and Kiihn 
120031 . 

A few laboratory measurements have been performed on 
amorphous solids in the FIR/mm range at low temperature. 
However, some of these studies were performed over a wide 



range of temperature and revealed a variation of the spectral 
behavior of the absorption coefficient with temperature. 

Agladze et al. ( 1996) performed absorption measurements 
on typical interstellar analog grains (crystalline and amorphous 
silicate grains) at low temperature ( 1 .2 - 30 K) in the millimeter 
range (700 /mi- 2.9 mm). They found an unusual behavior of 
the millimeter absorption of amorphous silicates : the absorp- 
tion coefficient decreases with temperature down to about 20 K 
and then increases again with decreasing temperature. They de- 
scribed the frequency dependence of the absorption coefficient 
using a temperature-dependent spectral index. Depending on 
the samples, the spectral index can decrease from f3 = 2.6 at 
10 K down to p — 1.8 near 30 K. They referred to tunneling 
effect in Two Levels Systems to explain this behavior. 

Mennella et al. d 19981 ) measured the temperature depen- 
dence of the absorption coefficients between 295 K and 24 K 
and for wavelength between 20 /mi and 2 mm on different 
kind of cosmic grain analogs (amorphous carbon grains, amor- 
phous and crystalline silicates ). They reported a significant 
temperature-dependence of the spectral index of the absorption 
coefficient, particularly strong for their amorphous iron-silicate 
sample. Their measurements showed a systematic decrease of 
the spectral index with increasing temperature. 

Boudet et al. (2005 ) performed measurements on different 
types of amorphous silicates (typical analog grains and SiC>2 
samples) for temperature between 10 K and 300 K and wave- 
length from 100 /mi to 2 mm. They found a strong temperature 
and frequency dependency of the absorption coefficient. They 
defined two spectral index depending on the wavelength range 
considered. For wavelengths between 500 //m and 1 mm they 
found a pronounced decrease of the spectral index with increas- 
ing temperature whereas, for wavelengths between 100 /mi and 
200 //m, the spectral index was found to be constant with tem- 
perature. They put some SiC>2 sample through a strong ther- 
mal treatment to remove most SiOH groups, and observed that 
the temperature-dependent absorption disappears totally. They 
thus identified the silanol groups as the defects that, in their 
silicate sample, are at the origin of the behavior. Considering 
that the OH groups can not simply increase the coupling be- 
tween the photon and the bulk phonons, they assumed reason- 
able that the defects induce closely spaced local energy min- 
ima, which correspond to the dynamics of a "particle" in an 
asymetric double-well potential. 

These laboratory studies provide strong indications that the 
absorption coefficients of amorphous grain analog materials 
vary substantially both with temperature and frequency. 

4. A model to explain temperature and spectral 
variation of dust optical properties 

4.1. Modeling disordered structure 

Two mechanisms have to be considered when dealing with 
amorphous materials, in order to take disorder into account. 

First, one should consider the acoustic oscillation excita- 
tion based upon the interaction of solids with electromagnetic 
waves due to Disordered Charge Distribution (See Schlomann 
I19641 ). This effect takes place not only in amorphous media, 
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but also in disordered crystals like mixed and polycrystals, 
and in some monocrystals with inverse spinel structure, for ex- 
ample. This absorption is temperature independent. The DCD 
model introduces a correlation length l c which quantifies the 
scale over which the atomic structure of the material realizes 
charge neutrality. The absorption spectrum of such a structure 
presents two asymptotic behaviors. Towards short wavelengths, 
the absorption is characterized by a spectral index /3 = 2 and in 
the longest wavelength range, the spectral index tends towards 
/3 = 4. The location of the spectral range in which the transition 
between the two regimes is directly related to the correlation 
length. The DCD model is described in Sect. 14.21 

The second mechanism describes the disorder at atomic 
scale, as a distribution of tunneling states, which leads to a 
temperature-dependent optical absorption, and enables to ex- 
plain the temperature dependence of the spectral index ob- 
served in laboratory experiments. This model was originally 
elaborated to explain the unusual thermal and optical proper- 
ties of amorphous material at low temperatures and has been 
described in details by Phillips ( 1972 ), Anderson, Halperin & 
Varma (119721 ) (see also Phillips 1987 for a review). In partic- 
ular, the TLS model was developed to explain the fact that, 
at low temperatures, the specific heat of amorphous solids ex- 
ceeds what is expected from the Debye theory for solids. This 
anomalous behavior is common to all amorphous materials and 
therefore appears independent of their detailed chemical com- 
position and structure. The existence of such two levels systems 
has been pointed out by means of microwave and submillime- 
ter spectroscopy experiments (e.g. Agladze et al. 119961 Bosch 
119781 and reference therein). The TLS model is described in 
Sect.gj] 

The two mechanisms considered here for the absorption by 
interstellar dust probably dominate in the FIR/mm and longer 
wavelengths domain. Both have a large degree of universality 
without specific signatures characterizing the chemical nature 
of the absorbing substance. On the other hand, they are sen- 
sitive to the internal structure of the solids, in particular their 
degree of ordering, and to mechanical properties such as den- 
sity and elasticity. 

4.2. DCD absorption 
4.2.1. DCD Theory 

In amorphous materials, the lack of long-range order permits 
single phonon/photon interactions with excitation of any modes 
of mechanical vibrations. Low frequency vibrations are not 
linked to molecular vibrational bands like is the case for ice 
or silicate bands. In infinite media they correspond to travel- 
ing acoustic waves. In finite bodies like interstellar dust grains, 
these modes can be described as standing waves. So we use the 
term "phonon" here for quantum of vibrational motion not re- 
stricting it (if not stated otherwise) to periodic lattice or infinite 
media. 

The concept of phonon quasiparticles (Tamm, 1930) arises 
as the result of quantization of vibrational motion in the har- 
monic approximation, which takes into account only quadratic 
terms (on atom displacements) of energy in Eq. d20l i. It corre- 



sponds to Eq. dT9b with bj = 0. Dissipation effects arise in the 
phonon concept as the result of phonon interaction, taking into 
account anharmonic terms in Eq. ( f2Qb . Their role was pointed 
first by Debye dl914| l. 

For harmonic oscillations u - ue~ ,a " with electrical field 
E(t) = Ee-' ,L " Eq. (O takes the form 



j Zj ;/ / m . 

re 1 



(34) 



where the coefficients of the dynamics matrix are 
j / j f 1 1 
The solution of the dispersion equation 



(35) 



- <./ /, - Dfj, /, = 

leads to resonant frequencies u> K and corresponding eigenvec- 

tors e K = (■■ ■ e , • • •)■ In an orthonormal basis of vectors e K , 

K j 

when 



KJ « J KK 

e j 

the equation d34"l > has a simple form 
{(j? K - a) 2 ) ti K — g K ■ E, 

where u K are amplitudes of normal oscillations, 



(36) 



(37) 



gi = Yj q i m i le 



and * denotes complex conjugates. 

The susceptibility corresponding to Eq. ( 1371 ) is given by the 
tensor 



r 1 . u)l - Co 2 - iury K ' 

KJ'J " 



(38) 



or in the macroscopically isotropic case by the scalar 
and therefore, 

1 aj 7Kqfq jmJ l e*fJ 
X = lm Xo = - y, — 2 JTTTT ^2' 



(39) 



Damping factors y K were introduced here to get a station- 
ary solution. y K values depends on phonon interactions, as dis- 
cussed above. In the FIR/mm region, overlapping of resonance 
curves produces a quasi continuous absorption spectrum. The 
integrated spectrum in Eq. d39b can be simplified replacing the 
integral of the sum of weak individual profiles by the sum of 
their integrated values. 



f 

xj — CX 



Im- 



n/2 



\w-ajf;\<Aa)/2 



-1 *£ f 
q jQ ,771 ■ e e ., 

1J 1J J KJ Kf 



(40) 



(41) 
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where Aa> is an averaging interval satifying 
a> » Ato » y K . 



(42) 



A number of states is proportional to Au>, therefore the re- 
sult does not depend on Aw and y K when condition in Eq. ( l42b 
is satisfied, which therefore also defines the validity region of 
the results of Vinogradov ( 1960) and Schlomann ( 1 1964b . 

In macroscopically uniform and isotropic cases the eigen- 
vectors of Eq. ( 1351 ) in the form of plane waves are 



(43) 



The manyfold of solutions of the dispersion equation is or- 
dered here in accordance with vector \k\ and splits into sev- 
eral branches labeled here by index r\. Each branch has a differ- 
ent dispersion curve u = a> n (\k\). Most branches are restricted 
to high frequencies (infrared) and only three branches extend 
down to the FIR/mm region. They are so called "acoustical" 
modes. For such low frequencies, these have a linear disper- 
sion dependance 



aj n = v n k, 



(44) 



where k - \k\, and v n is speed of sound. One branch corre- 
sponds to longitudinal propagation {r\ = /), the two others to 
2 polarizations of transverse propagation (rj = t). Numericaly, 
v, < v\. Displacements of neighbor ions in sound oscillations 
are synchronous and ej^ does not depend on j. The vectors 
of displacements = [e^, el, ei] are oriented such that 
£ki II k, ekt -L k. Their absolute values are defined by the nor- 
malization equation Eq. d36b . The number of normal oscillation 
per unit volume in the frequency interval Aa> is N^Aco, where 



N 0J 



\n 2 v\ 



and therefore 



cjN 
24^ 



i 2 II 

Vl \ m / oj/u, 

where N is number of ions per unit volume and 
If] - _L V 



(45) 



(46) 



■exp(i* •(/•;- /■;.))■ 



(47) 



The equation derived here is a generalization of the 
Schlomann's (1964) expression for one type of disordered ions 
in a perfect crystal (q 2 ) k m~ l , which was defined by a spatial 
power spectrum of charge distribution (q 2 )k only. 



4.2.2. Spectrum of DCD absorption 

The absorption spectrum corresponding to Eq. d47| > can be cal- 
culated using correlation functions. It is a cross-correlation 
function for general case 



c(Ar) = 



iqjqj 

\ 111; 



(48) 



Here ( ) designate an ensemble average, Ar = r,- - t? . 

Parameters q 2 v q 1 ^ are described below and do not depend on 
Ar. 

The correlation function can be determined from a physical 
model of considered material. The first, simplest model has an 
uncorrelated charge distribution and therefore a delta-function 
correlation 



4>(Ar) = 6(Ar) = 6(Ax)6(Ay)6(Az). 



(49) 



In general case a cross-correlation function, Eq. d49l gives 
a flat spatial spectrum with intensity not depending on k: 



1 , 

m Ik 

and therefore 



Xo(o>) 



12™? 



95 



(50) 



(51) 



Here, we omitted the second term in Eq. ( f46b since v] 3 <s 2u t 
An autocorrelation function should be used for m=const 



(52) 



< /> a (M = (q(r j )q(rd)(qi ) , 

with a flat spatial spectrum (q 2 )k = q\ for 6 correlation (Eq. 

leading to Eq. ( l5TT > with a constant q\ = q 2 /m. 

In perfect crystals the charge distribution is entirely corre- 
lated for long distances (long range order) and the absorption 
mechanism described by Eq. d46i i and Eq. (IBTt does not take 
place. Nevertheless even for crystals with a perfect lattice, a 
disordered charge distribution is possible if the lattice permits 
a non unique configuration of the charges in an elementary cell, 
as is the case, for example of the cubic lattice of salt (NaCl). 
In such cases, only stochastically distributed charges should be 
included in the calculation of q 2 , which can then be defined as 



7=1 



(53) 



where summing is to be done only for ions involved in the dis- 
ordered distribution and Aq denotes deviations from the mean 
charge value, which is different for different nodes of the crystal 
lattice. The ordered part of the charge distribution is linked only 
to the spectral components with a> > v n /d, where d is the lat- 
tice period. This frequency range corresponds to mid infrared 
absorption bands and are not considered here. 

In amorphous media the corresponding coefficient q 2 ^ 
should also exclude the regular part of the charge distribution 
produced by the short range order in the medium. An equation 
similar to Eq. ( |53l cannot be written in the general case. The 
effect of this short range screening can be expressed by the in- 
equality 

1 N a 2 

7=1 



The second model of Schlomann ( 1964) can be interpreted 
as screening over large scales. Large-scale electrical neutrality 
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of crystals and glasses could be a result of high charge mobility 
in solutions or melts from which the corresponding materials 
were manufactured. 

The corresponding correlation function should be written 
in generalized case as 

0(Ar) = 6(Ar) - -L exp(|Ar|// c ), (54) 

where l c defines the correlation length. The corresponding 
spectra are given by 

iq 2 )k = ~q\g(kl c ), (55) 
where 

g(x) = 1 - (1 + x 2 r 2 (56) 
and 

Xo(u) = T^— -T <7m g(.0)/0)c), (57) 

I2nv, 

where 



a 




Fig. 1. Absorption spectrum due to a Disordered Charge 
Distribution. The absorption coefficient q-dcd is given in arbi- 
trary units, normalized at A — 100/im and l c — lOnra. The 
solid, dot, dash, dash-dot and dash-dot-dot lines correspond to 
l c = 0.1, 0.3, 1, 3 and 10 nm respectively, assuming a transverse 
sound velocity v, — 3 • 10 5 cm/s. 



= v,/l c . 

The spectral dependence of Eq. 

u) 2 /a> 2 for to <s oj c 



is 



for to » io r 



(58) 



This is also true of other often used correlation function 
shapes, such as step correlation functions : 



<Kr) 



6(f) - 3/(4tt/3) for \r\ < l c , 







for \r\ > L 



(59) 



The spectral dependence in Eq. d58l is common for charge 
distributions, which are uncorrelated over short distances and 
respect neutrality over larger scales. 

The DCD absorption is temperature independent. The 
corresponding absorption coefficient a(u) ~ wx" (Eq.lT5]> 
presents two asymptotics behaviors on both sides of oj c which 
depends on the correlation length. In the high frequency range, 
the spectral index takes the value (3 — 2, and for the lower fre- 
quency range, the spectral index is equal to [5 — 4. This fre- 
quency dependence absorption is shown on Fig. [T] for various 
values of l c . 

4.3. TLS Absorption 

4.3.1. Phenomenological TLS Theory 

At low temperatures, the thermal and dielectric properties of 
disordered solids (e.g. glasses and mixed crystals) show defi- 
nite deviations from the predictions of the phonon theory devel- 
oped for perfectly ordered crystals. Most of these phenomena 
can be described in the formalism of the so called Two Level 
Systems (TLS), which is the simplest model of tunneling states. 
In comparison to DCD it could be considered as a second ap- 
proximation for the description of electromagnetic properties 



of solids linked to disordering. Vacancies unavoidable in dis- 
ordered lattices produce in first approximation a possibility for 
chaotic distribution of some charges, which was considered in 
the previous section. This distribution is considered as frozen 
in the lattice due to the large height of barrier dividing vacan- 
cies in comparison with the available thermal energy. The TLS 
theory takes into account transitions of charges due to quantum 
tunneling possible at any temperature. 

The TLS theory is the result of a macroscopic analysis of 
phenomena as opposed to an exact microscopic description. It 
is based on the hypothesis of a flat distribution of tunneling 
states with energy level differences sufficiently low with respect 
to the usual vibrational energy. The temperature dependance of 
two-level system properties results from deviations from the 
model of harmonic oscillator, which is a good approximation 
for the description of crystal lattice vibrations. The origin of 
TLS is vacancies in the lattice which gives the possibility, for 
atoms or groups of atoms, to have different equilibrium spatial 
configurations. As a result, the potential curve has at least a 
double-well form (DWP, see Fig. |2j in contrast to single-well 
potential (SWP) of harmonic oscillators in crystal without va- 
cancies. 

In the TLS model described e.g. by (Phillips 1987), the en- 
ergy distribution of states is defined under the assumption that 
the asymmetry A of DWP and the height V of the DWP poten- 
tial barrier are stochastic values and their probability distribu- 
tion can be assumed constant over a certain range of V and A 
values. Like for the well known tunneling probability, the en- 
ergy splitting Ao has an exponential dependence with barrier 
width d: 



A = HClexp[-d(2mV) 1/2 /H], 



(60) 



where SO = (E\ + Ej)j2 is the average energy of levels and 
m is the effective mass. This exponential dependence along 
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Table 1. Parameters of disordered solids. 



Material 


P 


£ 


u t 


k 


(q 2 )/e 2 




Pf4, 


7e 


Refs. 




g/cm 3 




cm/s 


nm 




D 




eV 




Silica Glass (Si0 2 ) 


2.2 


3.8 


4.1 • 10 5 


1-4 


1.1 


0.18" 


1 • 10" 5 


1.5 


1,2,3 


Soda-silica glass* (Na 2 O3Si0 2 ) 


2.9 


2.6 


3 • 10 5 


3 


1 




1.4 • 10- 3 


0.5 


4 


Chalogenic glasses 








1-4 


0.1-1 


0.5-5 






2 


Inverse spinel crystals 






4- 10 5 




1/14 








5 



a A previously reported value of dipole momentum 0.6 D was overestimated due to the absence of local field correction. 
b Barrier heights distribution parameters: maximum at V m /k B = 550 A", width V /k B = 410 A", V m j n /& B = 50 A". 
References: 1: Hubbard et al. l2TX)3l 2: Strom and Taylor [T9T7l 3: Phillips [19871 4: Bosch [T978l 5: Schlomann [T964l 




Fig. 2. Asymmetric double-well potential. Ordinate is the po- 
tential energy. Abscissa is a coordinate in configuration space, 
in the simple case a rotation angle or a linear displacement of 
groups of atoms in the lattice. A is the value of potential asym- 
metry. The difference in energies of the lowest levels results 
from the asymmetry A and tunnel splitting Aq. 

with the flat probability distributions of A and V, leads to the 
2-dimensional distribution 

P(A ,A)=?/Ao, (61) 

where P is a constant. The integral of P over Ao does not di- 
verge because the lowest Ao values are limited by the highest 
energy V through Eq.d60l>. 

The exact form of P(Ao) in the vicinity of A™ ln is not 
important in most cases because the corresponding energy 
(A™ n /fc B <K 0.01 K) are very far from values under consid- 
eration. A value of constant P is not denned by the TLS theory, 
it varies from one material to another, and can be regarded as a 
free parameter to be determined from experiments or observa- 
tions. 

In oder to better understand the TLS phenomenological ap- 
proach, Kiihn (2003) built a microscopic model using the trans- 
lational invariant Hamiltonian 

H = Z £ + Ui{Uj]) > 

j 



u = \ Z 3i ' (Ui - u i )2 + ^ Z (Mi " Uj)4 ' m 

•j 'j 

where pj are impulses, Uj are displacements of y'-th atom from 
some reference points. The first term is a quadratic potential 
with random coefficients Jy. The second term is a stabilizing 
quartic potential (not random) with a constant g > coefficient. 
The stochastically simulated coefficients 7 (J lead to disordered 
equilibrium arrangement of atoms. Disordered potential relief 
(Eq.l62l reproduces an ensemble of states both SWP and DWP, 
with a broad spectrum of barrier heights V and asymmetries A. 

Some TLS effects associated with excited-state transitions 
of TLS could be brought to the fore in some laboratory ex- 
periments (FitzGerald et al. 1994 and Sievers et al. 1998) on a 
number of mixed fluorite crystals and soda-lime silica glass. 
However these effects can only be seen in a narrow temper- 
ature range, roughly below 15K, before relaxation processes 
dominate at higher temperature. These effects will be incorpo- 
rated to our model in a forthcoming paper. In the following, we 
restrict our consideration to the two lowest levels with energy 

difference E = E 2 -E X = ^A 2 + Ajj. 

As opposed to thermal properties, the TLS absorption spec- 
trum depends not only on the TLS density of states distribution 
P over energies Ao, A, but also on the values of the matrix ele- 
ments of the dipole transitions between levels 

(l\p\2) = p h ^-, (63) 
E 

where p is an electric dipole moment, 
1 <3A 

«, = m (64) 

where £ is the value of the electrical field. 

The dipole moment p\, is the second TLS parameter, which 
characterizes its interaction with the electromagnetic field. For 
isotropic materials the square of dipole moment should be av- 
eraged over all directions as (p?) = p^/3. Like the parameter 
P discussed above, the parameter p\, varies from one material 
to another, and can be regarded as a free parameter, to be de- 
termined by experiments or observations. In some experimen- 
tal publications, values of the dipole moment are given, which 
are not corrected for the local field and orientation averaging, 
pt = p h (e' + 2)13 VS. 
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The TLS absorption spectrum shape can be obtained by 
solving the Bloch equations which describe the interaction of 
TLS with electromagnetic wave and lattice oscillations. The 
latter interaction is described by the third material dependent 
parameter of TLS theory, i.e. the elastic dipole moment given 
by: 



1 <9A 

7e= 2^' 



where e is the value of the strain field. The order of magnitude 
of fib and y e is about 1 D and 1 eV, respectively. 

The solution of general equations for the TLS absorption 
can be described as a result of two different contributions: an 
absorption which has a resonant character and an attenuation 
which has the typical form of a relaxation. In practice it is 
convenient to treat these processes separately. Bosch [T978] ex- 
plained his experimental results considering three processes: 
resonant tunneling , relaxation due to phonon-assisted tunnel- 
ing, relaxation due to phonon-assisted hoping over the poten- 
tial barrier. In the next sections, we will discuss each of them 
in more details. 

4.3.2. Resonant absorption 

The resonant absorption of a photon of energy tiu concerns 
only those TLS characterized by the splitting energy E2-E1 = 
ficj . Let N(a>) be the density per unit frequency of those TLS, 
and Ni{a>) and Nz(of) the densities per unit frequency of those 
TLS which are in the state of energy E\ and E2 respectively. Of 
course these densities are strongly temperature dependent, but 
in the TLS model only the two lower levels are considered and 
it is assumed that N - N\ + N2- 

The resonant tunneling absorption then takes the form 



The coefficient p g (Aq) S5 1 we introduced in Eq. (|68l l to de- 
scribe deviation of P(Ao) from Eq. (|6Tb at large Ao. To avoid 
confusion in this work we will not change definitions of vari- 
ables where it is possible, preferring transformation of the orig- 
inal equations taken from different sources. 

In the TLS theory, confirmed by low-temperature mea- 
surements of thermal properties, p c — 1 for energies smaller 
(65) than 10" 3 eV. If p Q = 1 for all Deltas, Eq. d68ll gives G(co) = 
PyU 2 =const. Expression for the TLS resonant tunneling absorp- 
tion, corresponding to p c = 1, was used by Bosch ( 1 19781 1 and 
Schickfus et al. fl975lll976l : 



tlLO 

c V? 



I e' + 2V 

7 B 4^— ) (M -#2). 



(66) 



where B\2 is the Einstein coefficient for absorption, 
47r 2 (l|/z|2) 2 



B 



3h 2 



The population of the levels under the thermal equilibrium 
hypothesis gives 

Ni - N 2 = N(oj)tanh(hoj/2k B T), 
and taking into account that E -tiu 
An 2 o) le' + 2\ 2 



3c V^ 7 



=(— ^— j tanh(/zw/2/fc B r) G(w), 



where 



G(cS) = 



Pp\ r tm p G (Ao)AodAo 
(M 2 Jo V 1 -(Ao/M 2 ' 



(67) 



(68) 



Equation d67|) was derived in the work of Hubbard et al. 
(2003), the term GicS) is the so called Optical Density Of States 
(ODOS) and can be determined from laboratory spectroscopy. 



4n 2 Pp' 2 

a res = =— to tanh(/jw/2A: B r). 

eye 7 



(69) 



The frequency independent ODOS was found to be 
in satisfactory agreement with experiments for wavenum- 
bers cj/2ttc <12cirT I . Its measured values {Pp\) are given 
in TableQ] Assuming p' = ID, the density of states is 
P = 0.9 ■ 10 33 erg -1 cirT 3 for soda-lime-silica glass. 

More accurate spectroscopic measurements in wide spec- 
tral regions show a definite drop of G(cl>) to zero within the error 
level, for u> higher than some frequency a> m . Agladze & Sievers 
( 1998 ) measured a profile with a cut-off at 20 crrT 1 in high den- 
sity amorphous phase ices of H2O and D2O. Fitzgerald, Sievers 
& Campbell (2001a) detected a sharp cut-off in fluoride mixed 
crystal spectra, with wavenumber <x> m l2nc =13 cirT 1 . Hubbard 
et al. (120031 1 observed a shallow cut-off at (x> m j2ixc =15 cirT 1 , 
9 cirT 1 and 6cirT 1 for the soda lime silica, the Si02 and the 
germanium glasses, respectively. 

Hubbard et al. (2003) suggested that the distribution 
P(A(>, A) (Eq.lrjTTi has a cut-off at energy E m = hu> m . In our 
Eq. ( f68b it corresponds to p c (Ao) in the form of a step function 



P G ( A o) 



1 for Ao < E n 
for Ao > E n 



(70) 



Under this assumption, Hubbard et al. (2003) calculated a 
frequency dependent expression for the ODOS by means of 
integration of distribution function equivalent to Eq. d68l l. For 
our consideration it is convenient to write it in the form 



G(w) = g(u/u) m )Pp h 



(71) 



The ODOS of (Hubbard et al. 2003) is described by the 
following function : 

[ 1 - -yl -(oj/a) m ) 1 forw/Wm > 1. 

However, the spectra corresponding to the ODOS proposed 
by Hubbard et al. (2003) (Eq.l72l exhibits a prominent break 
which was not observed in laboratory spectra of noncrystaline 
materials and has not been observed in spectra of astronomical 
dust emission. 

We suggest that the cut-off of the distribution function 
should be complemented by a continuity condition and there- 
fore continuity of p c (Ao). The simplest form of continuous 
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Fig. 3. Absorption spectrum due to resonant tunneling in Two 
Level Systems. The absorption coefficient is given in relative 
units normalized at A — lOOyum and T = \QK. The abscissa 
axis is normalized to the cut-off wavelength of the Optical 
Density of States A m = 700 pm (see text). The solid, dot, dash 
and dash-dot lines correspond to T — 10, 20, 50 and 100 K re- 
spectively. The various curves can be scaled using the simple 
relation /l m r=const. In the low frequency wings, curves scale 
as a ~ to 2 IT. 

function for p G (Ao) is considered in Appendix A and leads to 
the following a shape of the ODOS 



8(x) = 



1 + x 2 i 



for x < 1 , 



i (x) - gz(x) Vl - x^ 2 for x > 1 . 



(73) 



where x = co/co m , polynoms gi(x), gi(x) are defined in 
Appendix by Eq. ( |A.31 > and Eq. dA.4) . 

This ODOS function has no break at 0) = a> m , it is in better 
agreement with laboratory measurements of amorhous materi- 
als. It is also in full agreement with TLS theory at a> <K a> m . 
For this reason it will be used in the calculations of Sect. [5] 

4.4. Relaxation processes 

Direct consideration of the relaxation processes permits to eval- 
uate the nonresonant part of TLS absorption spectra through 
analysis of the mechanism of the TLS level population relax- 
ations after some deviation from equilibrium. The relaxation 
time constant r is determined by the rate (t _1 ) of transitions 
between levels, in which transitions due to intensive interac- 
tions with the lattice thermal oscillations (phonons) dominate. 

The rate of the TLS transitions as a result of strain field 
generated by thermal fluctuations in a lattice was evaluated by 
Phillips ([19721 . 



A 2 E- 2 , 



where 

T l =a £~ 3 tanh(E/2jfc B r) 
and 

5 i-4 -2 

a = nvjpn y e 



(74) 



(75) 



(76) 



is a material constant. Parameters included in Eq. d76*| i were al- 
ready described above. 

Phillips (1972) calculated also an equilibrium value of the 
TLS dipole moment Ji and therefore the resulting susceptibility 
X(cS) in the form of a Debye spectrum (Eq. 



#o'(«>) 



3k B TE 2 



sech 



i—)- 

\2JL771 



+ (wt) 2 



(77) 



The rate of tunneling relaxation (Eq.[74l> is proportional to 
the tunneling probability, so it does not take into account di- 
rect transitions in TLS which energy is higher than the bar- 
rier height. This transition rate is proportional to the number 
of phonons with such energy, which should follow a Boltzman 
distribution and depends on temperature and on the energy bar- 
rier height V. In glasses with V/k B k 500 K the time of hop- 
ping relaxation t varies with T < 100 K in accordance with the 
Arrhenius relation (Hunklinger& Schickfus 1981). 



t = T exp(y/^ B r) 



(78) 



with to ~ 10 -13 s (fused silica). 

Experimental values of V and To were determined using the 
temperature dependence of the position of the spectrum maxi- 
mum (a> = r _1 , see Eq. |77j. The distribution of TLS on V leads 
to the broadening of the corresponding spectrum, sometimes 
preventing definite detection of its maximum. The measure- 
ment of the spectrum broadening permits to evaluate the width 
and shape of the TLS distribution on V. In practice experimen- 
talists often use the shape of a truncated Gaussian probability 
distribution to fit the energy distribution of TLS. 



P(V) = C V Py(V), 



where 

Pv(V) 



i 



exp(-(V - V m ) 2 /V 2 ) 
\ 



for V > V min , 

for v < y min , 



(79) 



(80) 



where Cy is the usual normalization coefficient (Eq. ID.2b de- 
pending on V m , Vo, V m in, which are parameters of the distribu- 
tion shape. An example of measured numerical values of these 
parameters is given in TableQ] 

Starting from known distribution of TLS on energies 
E, Aq, V and spectrum (Eq.l77l> for fixed values of E, Ao, V 
one can synthesize a final relaxation spectrum. The integrated 
spectra are considered in the following sections. 

4.4.1. Tunneling relaxation 

The absorption spectrum due to phonon-assisted tunneling re- 
laxation can be obtained by integration of Eq. d77l i on E and Ao. 
The expression for this was provided by Fitzgerald, Campbell 
& Sievers.d2001bll 



ffphon - A - LO- F(C0, T), 

where 

3c Vi' 3 ) 



(81) 



(82) 
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Fig. 4. The function F p h on (p) (solid line) and the exponent of 
its slope P2 = iilnFphon/^ln P (dotted line) plotted as a func- 
tion of the parameter p = a p honw/(2fc B r) 3 (see text). The 
spectral index of phonon absorption y6 p hon is related to /?2 as 
P = 1 +/?2 and ranges from j8 P hon = 1 at low frequencies (low 
p values) to yS p h on = 2 at high frequencies (high p values). The 
temperature dependence p(T) generates a shift of the absorp- 
tion spectra along the frequency axis. 

is a material constant and : 



F(u,T) 2k K T X X, "l 1 " sech (^ B r)i 



u> dr dE 



+ C0 2 T 2 ' 

with T\ given by Eq. ( l75b . A change of variables made using 
relation between Ao and r for integral calculation. 

Another, simplified equation describing tunneling relax- 
ation absorption was used by Jackie (119721 1 and by Bosch 
< [[978) 



^phon 



f^sech^)- 
Jo 2k K T \2L.T/l 



at 2 T\ 



+ W 2 T? ' 



(83) 



where T\ is given by Eq. d75l >. parameters A and a (in Eq.|75]) 
are to be defined from experiments. The use of Eq. ( f83l is 
equivalent to evaluating the integral in Eq. d77b substituting val- 
ues by that for a symmetric DWP A=0. Formally it is equivalent 
to replacing VI - t\Jt by ti6(t - n). Our evaluation shows 
large inaccuracy in this procedure for values of absorption. For 
details see AppendixlBl 

However, this approach gives satisfactory accuracy for time 
constant evaluation linked to parameter a. For soda-silica glass 
(Bosch IT9781 I a = 4.2 • 10~ 56 erg 3 s. It corresponds to y e - 
0.5 eV (see Tabled). 

The full equation ( l83l for calculation of phonon relaxation 
spectra can be simplified without loss of accuracy. We made it 
using a function of one argument F v h on (p) instead of a function 
of two arguments F{a>, T) in Eq. (T8Tb . We introduce the inter- 
mediate parameter p{a>, T) = aw/(2A: B r) 3 . This permits to use 
a one-dimensional grid of argument, efficiently reducing com- 
puting time. The use of interpolation in precalculated values 




Fig. 5. Spectra of absorption due to tunneling relaxation a p hon 
in the FIR/mm. The absorption coefficient is given in arbitrary 
units, normalized at A = 100;um and T = 10 K. Line styles are 
the same as in Fig. [3] 



(Table lCTl in Appendix) make use of this function as fast as a 
standard one. 

Therefore an accurate and relatively simple equation for 
tunneling relaxation spectrum is 



= A ■ a) ■ F phon (a(L>(2k„T) ), 



^phon 

where A (Eq.| 
is given by 



(84) 



is a material constant and the function F p ho 



phon 



(P) 



/TV' 



tanx 



tanx 



dxdy, 



x\ - arctan(/? y arcth y) 



(85) 
(86) 



A plot of function F p hon(p) is shown on Fig. [4] 

The absorption spectrum for tunneling relaxation is tem- 
perature dependent in accordance with Eq. (fSTb . As opposed 
to resonant absorption, the tunneling relaxation absorption in- 
creases with temperature, as shown on Fig. [5] 

4.4.2. Hopping relaxation 

The expression for the absorption spectrum due to hopping re- 
laxation can be obtained by integrating Phillip's equation ( ITTl i 
over a Gaussian distribution of barrier height p y {V) given in 
Eq. (Ell : 



ffhop = Shop^) ^ 



1 + W 2 T 2 



P(V) dV, 



(87) 



where t = Toexp(V/k B T) and the coefficient Bho P (T) does not 
depend on u>. 

Examples of spectra described by Eq. d87l i are shown on 
Fig-H] It is known from general principles that the asymptotic 
behavior of e(w) when co — > leads to a(ai) ~ oj 2 (spectral in- 
dex y8=2). Figure|6]shows that this region begins at / < 10 MHz 
for T- 100 K and at /< 100 Hz for T < 50 K. 
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V (Hz) 

Fig. 6. Hopping relaxation absorption spectra in the frequency 
range from 1 Hz to 1 THz. The ratio of the absorption coef- 
ficient ahop to the parameter <yZ?hop/2 is plotted to show the 
whole spectrum at comparable scales and to evidence varia- 
tions of the peak position with temperature. The curves are 
computed using a barrier height of V m /fc B =550K (see Tab. [J}. 
Line styles are the same as in Fig. [3] Thin curves correspond to 
unbroadened spectra with Vo <K V m . Thick curves correspond 
to V /k B =41 OK. 



a/2-n (GHz) 




Fig. 7. Spectra of absorption due to hopping relaxation o-hop 
in the FIR/mm. The absorption coefficient is given in arbitrary 
units, normalized at A = 100 /mi and T = \QK. Line styles are 
the same as in Fig. [3] In this relatively narrow spectral region, 
the behavior of the curves is in agreement with the approxima- 
tion ofEq.©. 

At high frequencies (beginning from mm-wave region) the 
spectra shown on Fig. [7] are determined by the low energy part 
of the barrier height distribution. It can be approximated by 



COT 



!Vimn 1 + OJ 2 T 2 

where P m j n = P(Vmm), and therefore 



dV, 



• 

For spectra on Fig. [7] a value of P m j n corresponds to 
P min k s =exp[-(V mm - V m ) 2 /V 2 ]/V V^= 3.11-lO^K- 1 . 



(88) 



(89) 



Bosch ( 1 19781 1 used an approximated expression B hop = 
h r /T, where h r is a constant. The approximation used by 
Fitzgerald et al. (2001b) gives B hop = 2nfi 2 /3c yfe' which is 
a material constant. In Annex [Dj we derive a more accurate 
approximation which gives 



Bho P (T) 



in (e' +2) 2 - 2 



Pf4,(c A + hiT), 



3cV? 9 

where T is in Kelvin and material constant c A is given by 



(90) 
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(91) 



in accordance with Eq. d90T > and ( ID.4l i. c A is therefore an ad- 
ditional parameter of tunneling states. Its value corresponds 
to tunnel splitting in low A region where tunneling state den- 
sity distribution becomes significantly less than that of Eq. doTl i. 
The TLS theory was proved to apply in low temperature exper- 
iments at T down to 2 mK (Phillips [T9"8"7l >, therefore A™ in //t B < 
2 • 10~ 3 K and c A should verify c A ^ 5.8. The dependence of 
Bhop{T) is relatively weak and a tenfold variation in T from 
10 K to 100 K leads to variation of B hop (T) less than 30%. 

A difference in temperature dependence of considered ap- 
proximations is significant, especially for absorption calcula- 
tions over a wide range of temperatures. It does not manifest 
itself when fitting laboratory data over a limited region of tem- 
perature, leading only to a biased value of a parameter A. The 
main temperature dependence of a-hop is defined by exponential 
dependence of r on T (Eq.lTFIi. 

In the work of Bosch ( [19781 . Fitzgerald et al. d2Q01bl > and 
others, a comparison between experimental and calculated tem- 
perature dependancies was performed . It showed that relax- 
ation processes should be taken into account for temperatures 
T ^ 10 K. At these temperatures, relaxation dominates over 
resonant tunneling. Hopping relaxation becomes significant at 
T about 25 K and dominates at higher temperature. This ten- 
dency has a general character, and shows a similar behavior 
for substances as different as fluoride mixed crystals and silica 
glasses. 



5. Model of the dust FIR/mm emission 

The theoretical considerations given in Sect.|4]provide the for- 
malism for calculations of dust emissivity in the far-infrared 
and millimeter-wave region and can be gathered into a new 
model for explaining the astronomical emission in this spectral 
range. So far, modeling the astronomical dust emission spec- 
tra in this spectral range has been performed on purely phe- 
nomenological grounds with little connection to solid states 
physics. Unlike previous models, the new approach proposed 
here does not assume that dust is composed of crystalline ma- 
terial, but instead uses for the first time theoretical results ap- 
plicable to disordered materials which are likely to compose 
astronomical dust particles. 
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5.1. Model description 



The mass absorption coefficient of dust entering Eq.Q~|can be 
expressed as 



Kdcd = #m(3i/ f 3 cp) l , 



Kdust = K m + /fDCD + ^T L s- 



(92) 



The first term k 1r includes the opacity in infrared bands and 
other relatively short-wave absorption processes. References to 
dust models describing k 1r are given in Sect.Q] IR bands in the 
dust spectrum are linked to the optical modes of lattice oscilla- 
tions, like bending vibration bands of interatomic bonds such as 
O-Si-O or H-O-H. As was shown in the previous sections, the 
long wavelengths wing of infrared bands is relatively weak in 
the FIR/mm. Our model does not provide a new description of 
these processes and is complementary to IR models. While k 1r 
depends mainly on the chemical composition of dust, the fol- 
lowing terms kdcd and k jls are linked to the disordered struc- 
ture of the grain material and are largely independent of the 
dust chemical composition. 

As described in Sect.[4] kdcd and k tls are linked to the ex- 
istence of disorder in the structure of the grain material. Given 
the fast decrease of k 1r with wavelength, these effects are likely 
to dominate absorption in the FIR/mm. Disordered media ab- 
sorption takes place in all substances excluding perfect crystals. 
Amorphous dust as well as partially ordered or dirty crystalline 
dust can therefore be described by the model. 

The first term, kdcd corresponds to disordered charge ab- 
sorption described in Sect. l4.2l It is temperature independent. 
The second term, k jls is connected with transition of charges 
in a lattice between vacant potential minima due to tunneling 
or thermal activation (Sect. l4.3l ). This term displays a specific 
spectrum of absorption over a wide range of frequencies a>/2n, 
from about 1 Hz to about 1 THz. It includes three components 
designated in the model as resonant (k ks ), tunneling relaxation 
("■phonX and hopping relaxation (Vhop) respectively : 



< TLS ((0, T) = K res (w, T) + K phon (a>, T) + K hop (a>, T). (93) 



Spectral index of k tls significantly deviates from j3 — 2 and 
depends on temperature. Resonant opacity /c res decreases with 
T, relaxation opacity K p h on and Kh op rises with T. 

The various opacities of the model can be summarized as, 



KDCd(w) - KftCDtO 



1 - (l + a?lu?^) 

(oj \ i Ha) \ 
— Itanhl — I, 

fphon(w, T) - ^3^TLS W -fphon(flphonWr 3 ), 

i ts 2 v < ^ ts r ^(y,T)P(V)dV 
K bop (a>,T) = -K TLsW (c A+ lnD j __ 5 - 5 __,(97) 



(94) 
(95) 
(96) 



where t(V, T) = to exp(V/A: B r). 

The amplitude factors for the DCD and TLS terms can be 
expressed with respect to material properties as follows: 



7m 

4 ra 



K JLS = -n*P»l (cp)-\ 
flphon = ny^ 2 v 5 t ph A {2k^) 



(98) 
(99) 
(100) 



The function g KS is the normalized ODOS spectrum given 
by Eq. j73l . The function F re \ is defined by Eq. ( IC.lt and can 
be calculated interpolating Eq. ( |C.3t using precalculated values 
in Table lC.ll of the Appendix. 

The Gaussian probability distribution function of the TLS 
potential barrier height P(V) is defined by Eq. ( f80b with param- 
eters V m , V Q and V min . 

Kdcd and ^ TLS define the amplitude scale of kdcd 
and k tls . Parameters io c , io m , a p hon, t define the scale of 
*"dcd, f res , fphon, ^ho P , respectively, along the oj axis. 

The numerical values and units of the various parameters 
entering the model are given in Tab.[2] They set the frequency 
and temperature dependence of the model. 

It is important to note that the values proposed here only 
give an order of magnitude for the parameters, and are de- 
rived only from laboratory measurements (Sect. 13.231 and |4]l 
or theoretical expectations. They provide reference values only. 
They are not tuned to reproduce observed astronomical spectra. 
Derivation of the parameter values to be used for astronomical 
purposes will be derived in a forthcoming paper. 



Table 2. Parameter values. 



parameter 


value 


unit 




1.2- 1(T 24 


2 -1 

cm g s 


K + 

TLS 


1.5- 1(T 13 


cm 2 g _l s 


0)^ 


3.10 12 


s-' 




2.7 ■ 10 12 


s- 1 


^phon 


5.3 • 1(T 10 


K 3 s 


T ** 

' o 


io- 13 


s 


v m " 


7.6- 1(T 14 


erg 


Vo" 


5.7 • 1(T 14 


erg 


'min 


6.9 • 1(T 15 


erg 


C * 


5.8 





t values calculated using Eq. (I98H00) and the following approximate 



mean values 



3-10 cm/s, p = 3 g/cm 



p/4 



1CT\ y e = leV, 



= e 2 /m„, where e is the electron charge and m is the mass of the 
Oxygen atom. 

? using l c = v t /(o c =l nm (see Tab.flJ 

b corresponds to 2nc/ai m =0J mm (see Sect. l4.3.2l 

from TablefT] 

** T from Sect.l44l 

* evaluated from Sect. 14.4.21 



5.2. Model applicability 

5.2.1 . Wavelength range of validity 

The model presented here applies to FIR to millimeter emis- 
sion from amorphous materials. In the astronomical context, 
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it should be used to describe the emission of large dust grains 
which usually dominate the observed spectra in the FIR/mm. In 
principle, it should also apply to the emission from amorphous 
very small 3-D particles which emission usually dominates the 
mid-IR region of astronomical spectra. 

The DCD component is valid only in the frequency range 
dominated by acoustic oscillations. Towards high frequencies 
the corresponding limit is set where the sound wavelength is 
of the order of the interatomic distances d. For silicates this 
corresponds to an electromagnetic wavelength Ad lOjum. At 
shorter wavelengths IR models should be used. Towards the 
long wavelengths, the validity is limited to a> < co a with 

to a * TTUt/a, (101) 

where v, is the transverse sound velocity and a is the grain 
size. The DCD model should therefore be valid for all wave- 
lengths below A a — 3 mm for the range of sizes usually assumed 
for the BG component (e.g. 15-1 10 nm in Desert et al. 1990), 
using v, = 3 10 5 cm/s. Application of this part of the model to 
small transiently heated dust particles in the long wavelength 
range should therefore be handled with care. 

The TLS part of the model is in principle not limited to a 
given spectral range. 

5.2.2. Applicability in dust models 

In general, modeling FIR/mm dust emission from an interstel- 
lar source be seen as a three-step process. The first step is to 
select the materials which constitute the grains, and specify the 
optical properties (complex dielectric constants or refractive in- 
dex, or opacities) of the selected materials. The second step is 
to chose, for each dust component, a distribution of size, shape, 
composite structure, orientation with respect to the direction of 
the incoming radiation field, in order to calculate extinction, 
absorption and diffusion using Mie or some multipole approx- 
imation method. The third step is to specify the dust density 
and radiation field distributions, and to include radiative trans- 
fer to calculate the dust emission on a given line of sight. The 
model presented here concerns the first step of this dust emis- 
sion modeling only. 

The equations proposed in Sect.[5] do not include the lo- 
cal field correction factor given in Eq. (TT~5T >. They are valid for 
spherical or quasi-spherical particles. Any possible effect due 
to anisotropic shapes, such as a polarized emission by grains 
for example, is beyond the scope of this paper, and refers to the 
second step of modeling. 

The model equations do not depend on macroscopic values 
of the dielectric constant and average density. All the parame- 
ters involved in the model have a local, microscopic meaning. 
Therefore the derived emissivities are defined by the total mass 
of dust, and do not depend on the porosity of the grains. The 
model should therefore be equally valid for fluffy and for bulk 
grains. 

It is important to note that unlike standard modeling, our 
model shows that it is not sufficient for the first modeling 
step to select from databank some complex dielectric constants 
e(u>, T) obtained for bulk mineral samples of chemical compo- 
sition corresponding to interstellar element abundancies. First, 



the particle size should be taken into account. Indeed, if the 
opacity (mass emissivity) of dust (Eq.0 does not depend on 
particle size for a given dieletric constant e(w), the dieletric 
constante(w) itself is size dependent, as shown in Sect. 15. II 
Second, both DCD and TLS processes are more sensitive to 
the degree and the type of disorder in a given material, than to 
the exact chemical composition of the grain. 

5.3. Temperature and frequency dust opacity behavior 

In the general case, the model can produce a wide variety of 
spectral shapes, ranging from a "classical" spectrum /3-2 when 
only DCD is taken into account (Ktls = 0) and the correlation 
length caracterizing DCD is assumed infinite (o> c -Q) to spec- 
tra with much flatter behaviors when TLS effects are included, 
with mean spectral index values f3 over the whole FIR/mm 
range close to one at high dust temperatures. Some extreme 
cases even include the possibility of having /3 < or /3 > 2 
over some limited regions of the spectrum, due to the influence 
of the resonant TLS absorption and to the DCD cut-off respec- 
tively. 

However, the model generates a "typical emission profile" 
when a correlation length of nanometric scale is assumed, and 
when the TLS effects are expected and taken into account. This 
emission profile presents some interesting characteristic devi- 
ations compared to the commonly adopted profile (a modified 
black-body type emission characterized by a constant spectral 
index /3=2). To illustrate this fact, we show on Fig. [8] the dust 
opacity spectrum including both the DCD and TLS compo- 
nents (/Cdust = +kdcd + *r LS ), calculated using Eq. ( 19411971 ) and 
adopting the parameter values given in Tab.ff] First, the slope 
of the dust opacity versus wavelength, and consequently the 
spectral index, are no longer constant over the whole FIR/mm 
spectral range. As the temperature increases, it begins to de- 
viate in the longest wavelength range, starting from a local 
spectral index values higher than fi-2 around T=10K, towards 
lower values down to around /3=1. Second, as the temperature 
increases more, this spectral behavior propagates towards the 
short wavelengths of the FIR/mm range. It is remarkable that 
only with model parameter values referenced in the solid state 
physics literature, the temperature-dependent spectral behavior 
of the modeled dust opacity is in qualitative agreement both 
with laboratory data in that temperature range (Sect. 13. 231 and 
13331 

For each individual components of the model, the tempera- 
ture and frequency behaviors have been described in Sect.|4] 
The DCD opacity is temperature independent and shows a 
quadratic frequency behavior (kdcd ~ (Jj2 ) f° r high frequen- 
cies and drops dramatically (o> < o> c ) at lower frequencies. 
The transition between the two regimes is set by the value of 
the correlation frequency (w c ), or equivalently the correlation 
length of the charge distribution in the grain material. We point 
out that, although the transition between these two regimes 
is very smooth, this distinctive behavior opens the possibility 
to measure the correlation length from astronomical emission 
spectra. Note that the DCD effect is the only one with such a 
rapid fall-off of the absorption with wavelength, which seems 
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10° 10' 



Fig. 8. Total dust opacity spectra /cdcd + k tls ■ The spectra are 
given in arbitrary units, normalized at A = lOOji/m and T = 
10 K. The dash-dot-dot line shows the temperature independant 
DCD dust opacity alone computed for l c — 1 nm. Line styles for 
the temperature dependant total effect are the same as in Fig. [3] 

a / a m 

10° 10"' 




10° 10' 

Fig. 9. Temperature dependence of the resonant tunneling and 
phonon induced tunneling relaxation opacity (See Fig. [3] for 
comparison). The spectra are given in arbitrary units, normal- 
ized at A — 100//m and T = \QK. Line styles are the same as 
in Fig. [3] 

to be measured toward some astronomical sources which ex- 
hibit large ft values, including values higher than 2. 

Resonant tunneling opacity takes place around a> < co m and 
dominates at low temperatures since k ks ~ u> 2 /T (Fig. [3j. As 
was shown in Sect. l4.21 the wavelength A m = 27vc/a) m is in the 
spectral range from 500 fim to 1 mm for silicates. 

The tunneling relaxation spectrum is closely linked to the 
resonant tunneling. Their added effect is shown on Fig. [9] The 
amplitude of the resonant feature decreases with rising temper- 
ature, while the amplitude of the relaxation wings increases. 

The hopping relaxation behaves as a distinct spectral 
feature in the case of narrow barrier height distributions 
(Sect. |4.4l2l Fig.|6]l. The amplitude of the effect rises with fre- 
quency as /(hop ~ w 2 , until it reaches a maximum at a> = 
Tq 1 exp(-Vo/£ B r). At low temperature the location of the max- 



imum is shifted to MHz frequencies. It reaches the FIR/mm 
region at T - 100 K, a value which depends upon Vo. The hop- 
ping relaxation in the case of wide barrier height distribution, 
creates a very wide spectral distribution (Fig.|6]l. The behavior 
is then /Chop ~ TojP with a frequency-dependent slope value, 
usually ft < 1 (Eq.[12 Fig. 0. 

At a given temperature and wavelength, the relative effi- 
ciency of the TLS effects compared to the underlying DCD ef- 
fect is governed by the density and efficiency of the optically 
active TLS in the material, and so by the value of the parameter 
Pfit,. This value should be related to the fluffiness of the mate- 
rial and to the types and densities of defects in its structure. 

Generally, the resonant part of the TLS effects is expected 
to dominate over the other relaxation processes at tempera- 
ture below 10K, while the TLS hoping mechanism should fully 
dominate above 30K. 

5.4. Model relevance 

Recent astronomical observations of the interstellar dust emis- 
sion in the FIR/mm wavelength range show some unexpected 
features which can only be related with difficulty to previous 
knowledge about dust optical properties in that range (Sect. 13. ll 
All data have to be explained with a blackbody emission modi- 
fied by a frequency power law, whose exponent is constant over 
that FIR/mm range, with a temperature independent asymptotic 
value equal to two. Indeed, very little is known about the optical 
properties of the dust amorphous state in that range. Often stan- 
dard synthetic dust permitivities or opacities are used, which fit 
perfectly the near and mid-infrared observed features on the 
basis of a sum of Lorentz-type resonant profiles. But in the ab- 
sence of reliable data and knowledge, the FIR/mm behavior is 
often only modeled by the vanishing long-wavelength wings of 
those infrared resonances. In addition some laboratory data on 
various silicates reveal additional absorption features, strongly 
temperature dependent in the FIR/mm range. 

It is clear that only specific processes with characteristic 
frequencies located in the FIR/mm region are needed to ex- 
plain non purely asymptotic features in this spectral domain. 
Two phenomena satisfying these requirements are known in 
solid state physics. First, acoustic oscillations with acoustic 
wavelengths less than grain size have frequencies in the range 
3 • 10 ll) - 3 • 10 12 Hz (Sect. 15. lT >. These acoustic oscillations are 
"optically" active (interact with electromagnetic field) in dis- 
ordered media (Sect. 14. 2t . Second, tunneling states in the dis- 
ordered solids have frequencies A //z falling in the frequency 
range 5 ■ 10 5 - 5 • 10 11 Hz (Sect.gT2l). 

These two phenomena are linked to disordered material 
composing amorphous solids and dirty crystalline compounds. 
This is also the case of interstellar grains. They dominate the 
FIR/mm emission of the highly amorphous interstellar dust 
(Sect.Ol. 

The transverse sound oscillations and tunneling between 
the potential minima in a disorder lattice includes rotation of 
groups of atoms, limited and influenced by interactions with 
their neighboring atoms. In the classification of motion into 
electronic, vibrational and rotational states used in fundamen- 
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tal quantum mechanics, this model refers to rotational motions, 
when the former dust emissivity models used only asymptotic 
behavior of electronic and vibrational processes. In that sense, 
the model proposed here does not replace the former ones, but 
complements them by taking into account the physical phe- 
nomena efficient in the low frequency range. 

Both processes are sensitive to the degree and the type of 
disorder, rather than to the exact chemical composition of the 
material. As a consequence, the FIR/mm emission profile is 
governed by a few key-parameters which characterize the dis- 
order. There is no contradiction between this and laboratory re- 
sults of FIR/mm dust absorption. It is important to investigate 
the variations of the model parameters in laboratory measure- 
ments of various duly characterized samples. For this purpose, 
the model gives new avenues for such simulations, and such 
disorder characterization. It is likely that the FIR/mm proper- 
ties will be more influenced by the type and degree of disorder 
than by the chemical composition. 

It is important to note that the parameter values adopted in 
this paper were chosen from purely theoretical grounds or from 
laboratory results obtained on materials that may not be exactly 
representative of ISM dust particles. We make no attempt in 
this paper to model the actual astronomical observations. This 
will be performed in a companion paper. It is important that 
only the later parameters be used for astronomical purpose. As 
the type and degree of disorder in interstellar dust should be 
representative of the amorphous dust formation itself, we hope 
that this model will open the way to bring new insights about 
the interstellar medium dust. 

6. Conclusion 

The continuum FIR/mm emission is generally attributed to 
an interstellar dust component of amorphous grains, silicate- 
based, large enough 15-100 ram) to radiate at thermal equi- 
librium. So far the conception of FIR/mm dust emissivity was 
built on the semi-classical Lorentz or Drude models. They lead 
in the long wavelength limit to a dust emissivity obeying to 
a frequency power law, whose exponent (the spectral index) 
equals 2 and is usually considered as temperature independent. 
However a growing volume of FIR/mm observational data in- 
duces some doubts about the applicability of such models in 
that wavelength range. 

In this context we present a new model for FIR/mm dust 
emission, based on physical properties of the disordered mat- 
ter. We consider the interaction between the electromagnetic 
field and the acoustic oscillations in a Disordered Charge 
Distribution and a distribution of low energy Two Level 
Tunneling states. Both mechanisms apply to the amorphous 
materials with a large degree of universality, independant of 
the exact chemical nature of the dust. 

The proposed model predicts a disordered induced FIR/mm 
dust emission which dominates over the weak long wavelength 
wings of the infrared resonances. This emission is strongly 
temperature dependent, and leads to a millimeter emission en- 
hancement relative to predications from classical models. The 
FIR/mm emission spectra can not longer be characterized by a 
constant spectral index over that range. The detailed shape of 



the emission spectra is governed by a few key parameters which 
characterize disorder. The exact values of these parameters for 
interstellar dust should be determined from observations. 

Our model is complementary to the classical models which 
are reliable in the infrared, and adds emission processes due to 
disorder induced features. 
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Appendix A: TLS optical density of states 

The ODOS spectrum for continuous distribution of DOS is cal- 
culated by integration of Eq. (I(58l i over the energy plitting Ao. 
We use the following polynomial form of the distribution func- 
tion p c (x), with x = A()/fta> m 



. Ul + 3x 2 )(l -x 2 ) forx< 1, 
PoM=r o forx>l, 



(A.l) 



Equation (IA. lb is a least square polynomial fit to Eq. (l70l . 
removing the discontinuity at x=l. A similar result is obtained 
using the first Fourier harmonics of Eq. (ITOt . 

The calculated ODOS spectrum is given by G((l>) = 



Pjjlg(x) with x = aj/a) m and 



8(x) 



1 +x 2 gl (x) 



for x < 1 , 



1 + x 2 gi(x) - g 2 (x) Vl - x' 2 for x > 1, 



with 

g!(x)=4(5-6x 2 )/15, 

g 2 (x) = 8(l-x 2 )(2 + 3x 2 )/15.. 



(A.2) 



(A.3) 
(A.4) 



In the low frequency approximation (a> <k u> m ) we have 
p G ~ 1 and the ODOS function G(a>) does not depend on u>, in 
agreement with previous results of the TLS theory. 

Appendix B: Evaluation of simplified equation for 
tunneling relaxation calculation. 



A use of simplified single integral expression of Bosch (Eq.| 
instead of full double integral equation used here (Eq.[8TTl is 
based on assumption that income from asymmetric states is 
negligible. This approach gives correct result for resonant ab- 
sorption (Eq.l68l) where integrand term (1 - Arj/Sw) ' 2 am- 
plify weight of symmetric states Ao ~ ha>. Another situation 
takes place for relaxation absorption (Eq.[8Tb where integrand 
term (1 - ti/t) +1 ^ 2 emphasize weight of asymmetric states 
t » Ti . Direct calculations give discrepancy of results more 
than hundred times. A comparison could be made without de- 
tailed calculation be consideraton of noncoinsiding terms in in- 
tegrand expressions in integrals over dE. In case of (Eq.[83j it 

is f]_((Jj, t) = D 2 T\I(1 + Cl) 2 T 2 ). 

In case of (Eq.fSTb 



f 2 ((0,T)= Jl 

L V T (1 +OJ 2 T 2 



(B.l) 
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It is not difficult to define the asymptotics of both functions 



/l(W,Ti) 



U> T\ for U)T\ <K 1 
T7 1 for COT l » 1 



fl(OJ,Tl) 



{ am/2 for lot\ <c 1 
\ (2/3)t7' forwn » 1 

The asymptotic behaviors of f and f are different and the sim- 
plified equation is not accurate also for wide region of values 
of parameters. 

Appendix C: Calculation of integral for tunneling 
relaxation spectrum. 



The function F p hon(w, T) (Eq.[83ll can be simplified into a func- 
tion of a single parameter F p h on (p) with p = aa>/(2k B T) 3 
through a change of variables to the dimensionless variables 
x = arctan(dLtr) and y = tanh(£'/2fe B r). Equationl83lthen writes 

(ri-frX^s**; (Co 

Jo Jxi V tanx 



Fphon 

where 



x\ = arctan(/? y arcth 3 y), 



(C2) 



The function F^ on (p) can be further approximated using a 
fit of the form, 



(C3) 



Fphonip) = F phon (pi) \—\ , Pi < p < pt+i . 



Precalculated values of the function Fph on and the coeffi- 
cients /?2 entering the interpolated form are given in Tab. lC. 1 1 

Table CI. Values of F p hon(p;) an d Piipd entering the interpo- 
lated form of Fphon(p) in Eq. IC.3l 



Pi 


fphon (Pi) 


Piipd 


0.001 


1.4696 


-0.0321 


0.002 


1.4313 


-0.0446 


0.005 


1.3604 


-0.0678 


0.01 


1.2875 


-0.0924 


0.02 


1.1948 


-0.1248 


0.05 


1.0393 


-0.1826 


0.1 


0.8984 


-0.2397 


0.2 


0.7433 


-0.3090 


0.5 


0.5334 


-0.4185 


1 


0.3866 


-0.5114 


2 


0.2623 


-0.6084 


5 


0.1417 


-0.7343 


10 


0.0826 


-0.8208 


20 


0.0455 


-0.8950 


50 


0.0194 


-0.9659 


100 


0.0098 


-0.9897 


200 


0.0049 


-0.9969 


500 


0.0020 


-1.0000 


1000 


0.0010 





Appendix D: Calculations of integral for hopping 
relaxation 

Following the TLS formalism, the hopping relaxation spectrum 
can be computed as the integral of Eq. (|77| | over a distribution 
of TLS and barrier heights, taking into account that in contrast 
to tunneling relaxation, the time constant r depends on barrier 
height V according to Eq. d8Tb and does not depend directly on 
TLS parameters Aq and A: 



3k B T 
E 



Jo J-oo Jo A o 
A 2 ,/ E \ cjt 

^ sech (2XT)TT7^ dAodAdy - 



(D.l) 



For calculating this expression, the convention of 
Hunklinger and Schickfus (1981) is used, which allows to re- 
duce the triple integral to a single integration over distribution 
P(V)dV, where P(V) is taken to be Gaussian. This conven- 
tion suggests the independence of the distribution P(Aq, A) and 
P{V), which permits to separate the integrals and transforms 
Eq. ( ID.U into Eq. d87l i. in a form similar to that of Fitzgerald et 
al. (12001 bi t. In this case P(V) is an independent Gauss proba- 
bility distribution following the standard normalization 

I P(V) dV = 1 
Jo 

which defines the normalization coefficient Cy in Eq. d79i > 



C v = 



i_ E jv m -v^\ + i_ 

2 \ Vq 12 



(D.2) 



where the term in brackets is usually about unity. 
The coefficient Bhop(T) in Eq. d87b is equal to 



4tt (e' +2) 2 2- 2 



J Ar J Ar \2k B Tl\ E 2 A k B T 

4n (e' + 2) 2 2- , k T 

= — ~ — PplQn + Ci), 

C V7 9 3 A™ 11 

where the constant C\ is given by 



(D.3) 



C\ = In 4 — 1 + I In arcth x dx 
Jo 



-0.441. 



(D.4) 



The integral dD.3b over dAo has been determined precisely. 
Integrating over dE was performed assuming that terms of or- 
der A™ m /k B T and higher orders could be neglected, taking into 
account that A™ n //t B r <sc 1. 



C. Meny et al.: Far-infrared to millimeter astrophysical dust emission 



21 



References 

Agladze N.I. et al. 1996, ApJ, 462, 1026 

Agladze N.I., Sievers A.J. 1998, Phys. Rev. Lett., 80, 4209 

Agladze N.I., Sievers A.J., Jones S.A., Burlitch J.M., Beckwith S.V.W. 

1994, Nature, 372, 243 
Allamandola L.J., Tielens A., Barker J.R. 1985, ApJ, 290, L25 
Anderson P.W., Halperin B.I., Varma CM. 1972, Phil. Mag., 25, 1 
Andre P. et al. 1993, ApJ, 406, 122A 
Andriesse CD. 1974, A&A37, 257 
Barreiro, R. B. et al., 2006, MNRAS, in press 
Bennett, C. L. et al. 2003, ApJS, 148, 97 
Bernard J.-P. et al. 1999, A&A, 347, 640 
Bosch M.A. 1978, Phys. Rev., 40, 879 
Bosomworth D.R., 1967, Phys. Rev., 157, 709 

Boudet N. et al. 2002, in Chemistry as a Diagnostic of Star Formation, 
ed. C. Curry and M. Fich (NRC Press, Ottawa, Canada), 257 

Boudet N., Mutschke H., Nayral C, Jger C, Bernard J.-P, Henning 
T., Meny C.,2005, ApJ, 633, 272 

Bouchet, F. R., Gispert, R. 1999, New A, 4, 443 

Brucato J.R., G. Strazzulla et al. 2004, A&A, 413, 395 

Chandler C.J. et al. 1995, ApJ, 455L, 93 

Debye P. 1914, Ann. Physique 43, 49 

Desert F.-X., Boulanger F., Puget J.L. 1990, A&A, 237, 215 

Draine B.T., Anderson N. 1985, ApJ, 292, 494 

Draine B.T., Lee H.M. 1984, ApJ, 285, 89 

Draine B.T., Li A. 2001, ApJ, 551, 807 

Dupac X. et al. 2001, ApJ, 553, 604 

Dupac X. et al. 2002, A&A, 392, 691 

Dupac X., Bernard J.-P. et al. 2003, A&A, 404, Lll 

Dupac X., del Burgo C. et al. 2003, MNRAS, 344, 105 

Dwek E. et al. 1997, ApJ, 475, 565 

Ewald P. P., 1922, Naturwiss. 10, 1057 

Finkbeiner D.P., Davis M., Schlegel D.J. 1999, ApJ, 524, 867 

Fitzgerald S.A., Campbell J.A., Sievers A.J., 1994, Phys. Rev. Lett. 
73, 3105 

Fitzgerald S.A., Sievers A.J., Campbell J.A. 2001a, J. Phys. Condens. 
Matter, 13, 2177 

Fitzgerald S.A., Campbell J.A., Sievers A.J. 2001b, J. Phys. Condens. 
Matter, 13, 2095 

Galliano E, Madden S.C, Jones A.P., Wilson CD., Bernard J.-P. Le 

Peintre E, 2003, A&A, 407, 159 
Galliano, E, Madden, S. C, Jones, A. P., Wilson, C. D., & Bernard, 

J.-P. 2005, A&A, 434, 867 
Gordon M.A. 1988, ApJ, 331, 509 
Gordon M.A. 1990, ApJ, 352, 636 

Gromov V. et al. 2002, in Experimental Cosmology at mm-waves, ed. 
M. De Petris, M. Gervasi, Amer. Inst, of Phys. Conf. Proc. Ser., 
616, 205 

Hadni A., 1970, in S.S. Mitra and S. Nudelman (eds.), Plenum Press, 

New York, 561 
Henning T., Mutschke H. 1997, A&A, 327, 743 
Heuer A. 1997, Phys. Rev. Lett., 78, 4051 
Hubbard B.E. et al. 2003, Phys. Rev. B, 67, 144201 
Hunklinger S., Schickfus, M. v., 1981, in Amorphous Solids: Low 

Temperature properties, ed. WA. Phillips (Springer- Verlag), 81 
Jackie J. 1972, Z. Physik A, 257, 212 
Jager, CD. Fabian et al. 2003, A&A, 401, 57 
Kemper, E, W. J. Vriend, et al. 2004, ApJ, 609, 826 
Koike C, Hasegawa H., Manabe A. 1980. Ap&SS, 67, 495 
Kiihn R. 2003, Europhys. Lett., 62, 313 
Lamarre J.-M., 1994, Infrared Phys. and Techn., 35, 277 
Lamarre J.M. et al. 2003, New A Rev., 47, 1017 



Landau L.D., Lifshits E.M. 1982, Course of Theoretical Physics, vol. 

VIII (Moscow, Nauka) 
Leger A., Puget J.L. 1984, A&A, 137, L5 
Li A., Greenberg J.M. 1997, A&A, 323, 566 
Mathis J.S., Rumpl W„ Nordsieck K.H. 1977, ApJ, 217, 425 
Mennella V. et al. 1998, ApJ, 496, 1058 

Mon K.K., Chabal Y.J. and Sievers A.J., 1975, Phys. Rev. Lett., 35, 
1352 

Naselsky, P. D. et al„ 2005, Int. J. Mod. Phys. D, 14, 1273 
Oldham P. et al. 1994, A&A, 284, 5590 
Pajot E, 2006, A&A, 447, 769 
Phillips, W. 1972, J. Low Temp. Phys., 11, 757 
Phillips W. 1987, Rep. Prog. Phys., 50, 1657 
Rast H.E., Caspers H.H., Miller S.A. 1968, Phys. Rev. 171, 1051 
Reach W.T., Dwek E., Fixsen D.J., Hewagama T. et al., 1995, ApJ, 
451, 188 

Ristorcelli I. et al. 1998, ApJ, 496, 267 
Rubens H. and Hertz G., 1912, Berlin. Ber. 14, 256 
Schlomann E. 1964, Phys. Rev. A, 135, 413 
Schwartz P. 1982, ApJ, 252, 589 

Serra G. et al. 2002, Advances in Space Research 30, 1297 

Sievers A.J., Tu J., Agladze N., FitzGerald S.A., Campbell J.A., 1998, 

PhysicaB, 244, 159 
Stepnik B. et al. 2003 clim.conf., 1 87 Cross-Calibration of PRONAOS 

& ISO, in The Calibration Legacy of the ISO Mission (ESA SP- 

481) 

Stepnik B. et al. 2003, A&A, 398, 551 
Strom U., Taylor PC 1977, Phys. Rev., 16, 5512 
von Schickfus M., Hunklinger S., Pich L. 1975, Phys. Rev. Lett., 35, 
876 

von Schickfus M., Hunklinger S. 1976, J. Phys. C, Solid State Phys., 
9, L439 

Tamm I. E. 1930, Z. Physik, 60, 345 

Tegmark, M., Eisenstein, D. J., Hu, W, de Oliveira-Costa, A., 2000, 
ApJ, 530, 133 

Vinogradov, V. S. 1960, Fiz. Tverd. Tela, 2, 2622 (English transl. 

1961, Sov. Phys. - Solid St. 2, 2338) 
Walker C. et al. 1990, ApJ, 349, 515 
Weiland J.L. et al. 1986, ApJ, 306, L101 
Woody D. et al. 1989, ApJ, 337, 41 
Wright E. et al. 1992, ApJ, 396, L13 



